Grassmannian Frames with Applications to 
Coding and Communication 

Thomas Strohmer*and Robert Heath^ 



Abstract 

For a given class T of uniform frames of fixed redundancy we 
define a Grassmannian frame as one that minimizes the maximal cor- 
relation \{fk,fi)\ among all frames {fk}kei S T. We first analyze 
finite-dimensional Grassmannian frames. Using links to packings in 
Grassmannian spaces and antipodal spherical codes we derive bounds 
on the minimal achievable correlation for Grassmannian frames. These 
bounds yield a simple condition under which Grassmannian frames co- 
incide with uniform tight frames. We exploit connections to graph the- 
ory, equiangular line sets, and coding theory in order to derive explicit 
constructions of Grassmannian frames. Our findings extend recent re- 
sults on uniform tight frames. We then introduce infinite-dimensional 
Grassmannian frames and analyze their connection to uniform tight 
frames for frames which are generated by group-like unitary systems. 
We derive an example of a Grassmannian Gabor frame by using con- 
nections to sphere packing theory. Finally we discuss the application 
of Grassmannian frames to wireless communication and to multiple 
description coding. 
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1 Introduction 



Orthonormal bases are an ubiquitous and eminently powerful tool that per- 
vades all areas of mathematics. Sometimes however we find ourselves in a 
situation where a representation of a function or an operator by an over- 
complete spanning system is preferable over the use of an orthonormal basis. 
One reason for this may be that an orthonormal basis with the desired prop- 
erties does not exist. A classical example occurs in Gabor analysis, where 
the Balian-Low theorem tells us that orthonormal Gabor bases with good 
time-frequency localization cannot exist, while it is not difficult to find over- 
complete Gabor systems with excellent time-frequency localization. Another 
important reason is the deliberate introduction of redundancy for the purpose 
of error correction in coding theory. 

When dealing with overcomplete spanning systems one is naturally lead 
to the concept of frames [TU]. Recall that a sequence of functions {fk}k£i 
is a countable index set) belonging to a separable Hilbert space TL is said to 
be a frame for TL if there exist positive constants (frame bounds) A and B 
such that 

kei 

for every / ETC. 

Even when there are good reasons to trade orthonormal bases for frames 
we still want to preserve as many properties of orthonormal bases as possible. 
There are many equivalent conditions to define an orthonormal basis {ek}k&i 
for TC, such as 

/ = X)</> e ^ V ^ e H > and H e *H = !> I, (2) 

kex 

or 

{efc}fc 6 j is complete in TC and (e k , e t ) = 8 k j, (3) 

where Sk,i denotes the Kronecker delta. 

These two definitions suggest two ways to construct frames that are "as 
close as possible" to orthonormal bases. Focusing on condition (J2J) we are 
naturally lead to uniform tight frames, which satisfy 

f = j y E,(fJk)h i V/€«, and||/ fc || 2 = l, Vkel, (4) 
kex 
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where A is the lower frame bound. This class of frames has been frequently 
studied and is fairly well understood [TT H lo] 12%) I2T] . 

As an alternative, as proposed in this paper, we focus on condition (jSJ), 
which essentially states that the elements of an orthonormal basis are per- 
fectly uncorrelated. This suggests to search for frames {fk}kei such that 
the maximal correlation \(fk, fi)\ for all fc,I 6 I with k ^ I, is as small as 
possible. This idea will lead us to so-called Grassmannian frames, which are 
characterized by the property that the frame elements have minimal cross- 
correlation among a given class of frames. The name "Grassmannian frames" 
is motivated by the fact that in finite dimensions Grassmannian frames co- 
incide with optimal packings in certain Grassmannian spaces as we will see 
in Section |21 

Recent literature on finite-dimensional frames fZT\ EI] indicates that 
the connection between finite frames and areas such as spherical codes, alge- 
braic geometry, graph theory, and sphere packings is not well known in the 
"frame community". This has led to a number of rediscoveries of classical 
constructions and duplicate results. The concept of Grassmannian frames 
will allow us to make many of these connections transparent. 

The paper is organized as follows. In the remainder of this section we 
introduce some notation used throughout the paper. In Section |21 we focus 
on finite Grassmannian frames. By utilizing a link to spherical codes and 
algebraic geometry we derive lower bounds on the minimal achievable corre- 
lation between frame elements depending on the redundancy of the frame. 
We further show that optimal finite Grassmannian frames which achieve this 
bound are also tight and certain uniform tight frames are also Grassmannian 
frames. We discuss related concepts arising in graph theory, algebraic ge- 
ometry and coding theory and provide explicit constructions of finite Grass- 
mannian frames. In Section |3] we extend the concept of Grassmannian frames 
to infinite-dimensional Hilbert spaces and analyze the connection to uniform 
tight frames. We give an example of a Grassmannian frame arising in Gabor 
analysis. Finally, in Section 0] we discuss applications in wireless communi- 
cation and coding theory. 

1.1 Notation 

We introduce some notation and definitions used throughout the paper. Let 
{fk}kex be a frame for a finite- or infinite-dimensional Hilbert space 7i. Here 
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X is an index set such as Z, N or {0, . . . ,N — 1}. The frame operator S 
associated with the frame {fk}kei is defined by 

sy = X)</, /*>/*• ( 5 ) 

k£l 

S is a positive definite, invertible operator that satisfies AI < S < BI, where 

1 is the identity operator on 7i. The frame analysis operator T : 7i — > ii^X) 
is given by 

Tf = {(f,f k )} k ^, (6) 
and the frame synthesis operator is 

T* :i 2 {X)^H: T{c k } kex = ^ c k f k . (7) 

feex 

Any f EH can be expressed as 

feex kei 

where {h k } k£ j is the canonical dual frame given by h k = S^ 1 f k . If A — B 
the frame is called tight, in which case 5* = A I and h k = \fk- The tight 
frame canonically associated to {f k } k ei is f k . 

If 1 1 /fell = 1 for ah k then {f k } k& j is called a uniform frame. Here 
|.|| denotes the ^-norm of a vector in the corresponding finite- or infinite- 
dimensional Hilbert space. Uniform tight frames have many nice properties 
which make them an important tool in theory [H3J [2Ej and in a variety of 
applications HOI EHl Ej . Observe that if {f k } k( =i is a uniform frame, then 
{S~z fk}kex is a tight frame, but in general no longer uniform! 

We call a uniform frame {f k }k& equiangular if 

|(/ fc ,/,)|=c for all fc, Z with k ^ I, (9) 
for some constant c > 0. Obviously any orthonormal basis is equiangular. 

2 Finite Grassmannian frames, spherical codes, 
and equiangular lines 

In this section we concentrate of frames {f k } k=l for E m where E = 1 or C. 
As mentioned in the introduction we want to construct frames {f k } k= i such 
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that the maximal correlation \{fk, fi)\ for all fc,l£l with k ^ I, is as small 
as possible. If we do not impose any other conditions on the frame we can 
set N = m and take {fk}k=i ^° an orthonormal basis. But if we want to 
go beyond this trivial case and assume that the frame is indeed overcomplete 
then the correlation \{fk, fi)\ will strongly depend on the redundancy of the 
frame, which can be thought of as a "measure of overcompleteness" . Clearly, 
the smaller the redundancy the smaller we expect \(fk, fi)\ to be. In E m the 
redundancy p of a frame {fk}k=i is defined by p = — . 

Definition 2.1 For a given uniform frame {fk}k=i in E m we define the max- 
imal frame correlation M.{{fk]k=i) by 

M({f k }ti)= m«{KA,/,>|}. (10) 

The restriction to uniform frames in the definition above is just for con- 
venience, alternatively we could consider general frames and normalize the 
inner product in (fTUj) by the norm of the frame elements. Hence without 
loss of generality we can assume throughout this section that all frames are 
uniform. 

Definition 2.2 A sequence of vectors {uk} k= i in E m is called a Grassman- 
nian frame if it is the solution to 

min{M({f k }ti)}, (11) 
where the minimum is taken over all uniform frames {fk}k=i ^ n ^ m ' ■ 

In other words a Grassmannian frame minimizes the maximal correlation 
between frame elements among all uniform frames which have the same re- 
dundancy. Obviously the minimum in (jllj) depends only on the parameters 
iV and m. 

Two problems arise naturally when studying finite Grassmannian frames: 
Problem 1: Can we derive bounds on M. ({fk}k=i) f° r given iV and ml 
Problem 2: How can we construct Grassmannian frames? 

The following theorem provides an exhaustive answer to problem 1. The 
theorem is new in frame theory but actually it only unifies and summarizes 
results from various quite different research 
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Theorem 2.3 Let {fk})? = i be a frame for W 71 . Then 



Equality holds in f!12j1 z/ and onfo/ i/ 

{/fcjfcLi is an equiangular tight frame. (13) 

Furthermore, 

Tfl\Yfl -I- 1 ) 

ifE = R equality in can only hold if N < — >— ^ (14) 

i/E = C tnen equality in (|12j) can on/y hold if N < m 2 . (15) 

Proof: A proof of the bound (*T2"|) can be found in [4*3*1 . It also follows 
from Lemma 6.1 in ["IT]. One way to derive (|12|) is to consider the non- 
zero eigenvalues Ai, . . . , A m of the Gram matrix R = {(fk, fi)}ki=v These 
eigenvalues satisfy YlT=i ^k = N and also 

m N N 2 

E^EEKA./.)!"^. (16) 

k=l k=l 1=1 

see |35| I41|. The bound follows now by taking the maximum over all \ (fk,fi)\ 
in (fi*o"|) and observing that there are N(N — l)/2 different pairs (fk, fi) for 
k, I G X with k^l. 

Equality in (|12j) implies Ai = • • • = A m = — , which in turn implies 
tightness of the frame, and also \(fk,fi)\ 2 = J^^-i) ^ or a ^ ^ wr th k ^ I 
which yields the equiangularity (cf. also [HI3"!]). Finally the bounds on N 
in (JI1J), (HU) follow from the bounds in Table II of |TT]. □ 

We call uniform frames that achieve the bound (fT2*|) optimal Grassman- 
nian frames. The following corollary will be instrumental in the construction 
of a variety of optimal Grassmannian frames. 

Corollary 2.4 Let m, N e N with N > m. Assume R is a hermitian N x N 
matrix with entries & = 1 and 
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for k, I = 1, . . . , N; k ^ I. If the eigenvalues Ai, . . . , Ajv of R are such that 
Ai = • • • = A m = — and A m+ i = • • • = Ajv = 0, then there exists a frame 
{fk}k=i i> n ^ m that achieves the bound (fT2*j) . 

Proof: Since R is hermitian it has a spectral factorization of the form R = 
WAW*, where the columns of W are the eigenvectors and the diagonal matrix 
A contains the eigenvalues of R. Without loss of generality we can assume 
that the non-zero eigenvalues of R are contained in the first m diagonal 

entries of A. Set fa := y~^{Wk,i}%Li for k — 1, . . . , N. By construction 

we have (f k ,fi) = Ri,k, hence {fk}k=i is equiangular. Obviously {fk}k=i is 
tight, since all non-zero eigenvalues of R are identical. Hence by Theorem 12. HI 
{fk\k=i achieves the bound (|T2|). □ 
On the first glance Corollary 12.41 does not seem to make the problem of 
constructing optimal Grassmannian frames much easier. However by using 
a link to graph theory and spherical designs we will be able to derive many 
explicit constructions of matrices having the properties outlined in Corol- 
lary o 

While the concept of Grassmannian frames is new in frame theory there 
are a number of related concepts in other areas of mathematics. Thus it is 
time to take a quick journey through these areas which will take us from 
Grassmannian spaces to spherical designs to coding theory. 

Packings in Grassmannian spaces: 

The Grassmannian space Q (m, n) is the set of all n-dimensional subspaces 
of the space W 71 (usually the Grassmannian space is defined for R only, al- 
though many problems can be analogously formulated for the complex space). 
Q{m, n) is a homogeneous space isomorphic to 0(m) / (0(n) x Oim — n)), it 
forms a compact Riemannian manifold of dimension n(m — n). 

The Grassmannian packing problem is the problem of finding the best 
packing of n-dimensional subspaces in E m , such that the angle between 
any two of these subspaces becomes as large as possible [7||H]. In other words, 
we want to find N points in Q (m, n) so that the minimal distance between 
any two of them is as large as possible. For our purposes we can concentrate 
on the case n — 1. Thus the subspaces are (real or complex) lines through 
the origin in E m and the goal is to arrange Af lines such that the angle be- 
tween any two of the lines becomes as large as possible. Since maximizing 
the angle between lines is equivalent to minimizing the modulus of the inner 
product of the unit vectors generating these lines, it is obvious that find- 
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ing optimal packings in Q{m, 1) is equivalent to finding finite Grassmannian 
frames (which also motivated the name for this class of frames). 

By embedding the Grassmannian space Q(m,n) into a sphere of radius 
y^n(m — n)/m in M. d with d = (m + l)m/2 — 1, Conway, Hardin, and Sloane 
are able to apply bounds from spherical codes due to Rankin to de- 
rive bounds on the maximal angle between N subspaces in Q(n,m), see the 
very inspiring paper [7j. For the case n — 1 the bound coincides of course 
with (JT2J). 

Spherical codes: 

A spherical code S(m, N, s) is a set of N points (code words) on the Tri- 
dimensional unit sphere Q m , such that the inner product between any two 
code words is smaller than s, cf. jS]. By placing the points on the sphere as 
far as possible from each other one attempts to minimize the risk of decoding 
errors. Antipodal spherical codes are spherical codes which contain with each 
code word w also the code word —w. Clearly, the construction of antipodal 
spherical codes whose N points are as from each other as possible is closely 
related to constructing Grassmannian frames. 

In coding theory the inequality at the right-hand side of is known as 
Welch bound, cf. Therefore uniform tight frames are known as Welch 

bound equality (WBE) sequences in coding theory 1 . WBE sequences have 
gained new popularity in connection with the construction of spreading se- 
quences for Code-Division Multiple-Access (CDMA) systems [121 I2H EE]- 
WBE sequences that meet f)12jl with equality are called maximum WBE 
(MWBE) sequences [HH EH]- While Welch (among other authors) derived 
the bound (JT2J) he did not give an explicit construction of MWBE sequences. 

Spherical designs: 

A spherical t- design 2 is a finite subset X of the unit sphere Q m in M m , 
such that 

\X\- lS ^h{x) = I h{x)dw{x), (18) 

for all homogeneous polynomials h G Hom t (lR m ) of total degree i in m vari- 
ables, see e.g. [HE]- A spherical design measures certain regularity properties 
of sets X on the unit sphere Q m . Another way to define a spherical t-design 
is by requiring that, for k — 0, . . . , t the k-th moments of X are constant with 

1 The authors of 021 incorrectly call WBE sequences tight frames. 

2 A spherical t-design should not be confused with an "ordinary" i-design. 
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respect to orthogonal transformations of M. m . Here are a few characteriza- 
tions of spherical t-designs that make the connection to the aforementioned 
areas transparent. For details about the following examples we refer to [T2"] . 
Let the cardinality of X be N. X is a spherical 1-design if and only if the 
Gram matrix R(X) of the vectors of X has vanishing row sums. X is a 
spherical 2-design if it is a spherical 1-design and the Gram matrix R(X) 
has only two different eigenvalues, namely N/m with multiplicity m, and 
with multiplicity N — m. An antipodal spherical code on Q m is a 3-design 
if and only if the Gram matrix of the corresponding set of vectors has two 
eigenvalues. 

Equiangular line sets and equilateral point sets: 

In |3U1 ITT] Seidel et al. consider sets of lines in M m and in C m having a 
prescribed number of angles. They derive upper bounds on the number of 
lines in the case of one, two, and three prescribed angles (in the latter case, 
one of the angles is assumed to be zero). Most interesting are those line sets 
that actually meet the upper bound. In 41 van Lint and Seidel consider 
a similar problem in elliptic geometry. Since the unit sphere in R m serves 
as model for the m — 1-dimensional elliptic space E m_1 where any elliptic 
point is represented by a pair of antipodal points in W 71 , the construction of 
equilateral point sets in elliptic geometry is of course equivalent to the con- 
struction of equiangular lines sets in Euclidean geometry. Recall that optimal 
Grassmannian frames are equiangular, hence the search for equiangular line 
sets is closely related to the search for optimal Grassmannian frames. 

Characterization of strongly regular graphs: 

Graphs with a lot of structure and symmetry play a central role in graph 
theory. Different kinds of matrices are used to represent a graph, such as the 
Laplace matrix or adjacency matrices What structural properties can 
be derived from the eigenvalues depends on the specific matrix that is used. 
The Seidel adjacency matrix A of a graph T is given by 

if the vertices x, y G T are adjacent, 
if the vertices x, y e T are non-adjacent, (19) 
if x = y. 

If A has only very few different eigenvalues then the graph is (strongly) 
regular, cf. [3]. The connection to Grassmannian frames {fk}k=i that achieve 
the bound (J 12)) is as follows. Assume that the associated Gram matrix R = 




9 



{(fk, fi)}ki=i has entries ±a and 1 at the diagonal. Then 

A=-(R-I) (20) 
a 

is the adjacency matrix of a regular two-graph [37] . We will make use of this 
relation in the next section. 

2.1 Construction of optimal Grassmannian frames 

In this section we give some explicit constructions for infinite families of op- 
timal finite Grassmannian frames. Note that optimal Grassmannian frames 
do not exist for all choices of m and N (assuming of course that N does not 
exceed (m + l)m/2 or m 2 , respectively). For instance there are no 5 vectors 
in M. 3 with maximal correlation In fact, although the 5 vectors in R 3 that 

minimize (jllj) are equiangular, the maximal inner product is (but not ^=), 
see |7]. On the other hand the 7 vectors in M. 3 that minimize yield a 
uniform tight frame, but not an equiangular one (which should not come as 
a surprise since the choice N = 7 exceeds the bound N < m(m + l)/2. Note 
that for C 3 we can indeed construct 7 lines that achieve the bound (|12|). We 
refer to [7j for details about some of these and other examples. 

Corollary 2.5 Let E = R or C and N — m + 1. {fk}k=i ^ s an optimal 
Grassmannian frame for E m if and only if it is a uniform tight frame. 

Proof: An optimal Grassmannian frame {fk\k=i with N = m + 1 can 
be easily constructed by taking the vectors to be the vertices of a regular 
simplex in E m , cf j7j. Thus by Theorem 12. 31 { fh}f 1 is a uniform tight frame. 
On the other hand it was shown in j2Ij that all uniform tight frames with 
N = m + 1 are equivalent. Since this equivalence relation preserves inner 
products it follows that any uniform tight frame {fk}k=i with N = m + 1 
achieves the bound (fT2|) . □ 

A uniform tight frame {fk}k=i with N = m + 1 also provides a spherical 
1-design, which can be seen as follows. When iV = m + 1 we can always 
multiply the elements of {fk}k=i by ±1 such that the Gram matrix R has 1 
as its main diagonal entries and — — else. Hence the row sums of R vanish 
and therefore {fk}k=i constitutes a spherical 1-design. It is obvious that the 
Seidel adjacency matrix A of a graph which is constructed from a regular 
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simplex has A^k = and A^i = — 1 for k ^ I, which illustrates nicely the 
relationship between A and R as stated in (|T9j). 

The following construction has been proposed in [7] . An nxn confer- 
ence matrix C has zeros along its main diagonal and ±1 as its other entries, 
and satisfies CC T = (n — l)I n , see Conference matrices play an im- 

portant role in graph theory j37j. If C*2 m is a symmetric conference matrix, 
then there exist exist 2m vectors in W 71 such that the bound (|12j) holds with 
p(2m,m) = 1/ ' \j2m — 1. If C^m, is a skew-symmetric conference matrix (i.e., 
C = —C T ), then there exist exist 2m vectors in C m such that the bound ()12|) 
holds with p(2m,m) = \ j \/2m — 1, see Example 5.8 in [TI]. The link be- 
tween the existence of a (real or complex) optimal Grassmannian frame and 
the existence of a corresponding conference matrix Ci m can be easily as seen 
as follows. Assume that {fk}k=i achieves ()12|) and denote a := l/y/2m — 1. 
We first consider the case E = R. Clearly the entries of the 2m x 2m Gram 
matrix R = {(fk, fl)}ki=i are -^M = i° ior k ^ I and i?*^ = 1. Hence 

C:=-{R-I) (21) 

a 

is a symmetric conference matrix. For E = C we assume that Rk,i = ±iu for 
k ^ / and Rk,k = 1- Then 

C:=— (i?-/) (22) 

is a skew-symmetric conference matrix. 

The derivations above lead to the following 

Corollary 2.6 (a) Let N = 2m, with N = p a + l where p is a prime number 
and a EN. Then there exists an optimal Grassmannian frame in M. m which 
can be constructed explicitly. 

(b) Let N = 2m, with m = 2 a with a G N. Then there exists an optimal 
Grassmannian frame in C m which can be constructed explicitly. 

Proof: Paley has shown that if iV = p a +l with p and a as stated above, then 
there exists a symmetric N x N conference matrix, moreover this matrix can 
be constructed explicitly, see (221 CHI- For the case N = 2m = 2 a+1 a skew- 
symmetric conference matrix can be constructed by the following recursion: 
Initialize 

(23) 



-1 

1 
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and compute recursively 



a 



2 m 



C m 



a — t 

^m ± n 



(24) 



then it is easy to see that C*2 m is a skew-symmetric conference matrix. 

An application of Corollary 12.41 to both, the symmetric and the skew- 
symmetric conference matrix respectively, completes the proof. Q 

Hence for instance there exist 50 equiangular lines in IR 25 with angle 
acos(l/v^49) and 128 equiangular lines in C 64 with angle acos(l/Vl~27). The 
construction in ()23jl . (|24jl is reminiscent of the construction of Hadamard 
matrices. Indeed, C m + J m is a skew-symmetric Hadamard matrix. 



2.1.1 Nearly optimal Grassmannian frames 



Theorem 12.31 gives an upper bound on the cardinality of optimal Grassman- 
nian frames. If the redundancy of a frame is too large then it cannot achieve 
equality in (|12p. But it is possible to design Grassmannian frames whose 
cardinality slightly exceeds the bounds in Theorem 12. 3\ while their maximal 
correlation is close to the optimal value. For instance there exist frames 
{fk}k=i m C m where N = m 2 + 1, with maximal correlation M. = hU. In 
fact, these nearly optimal Grassmannian frames are unions of orthonormal 
bases, and the modulus of the inner products between frame elements takes 
on only the values and hL. We refer to |U EI] for details about these 
amazing constructions, which find an important application in the design of 
spreading sequences for CDMA j2*4~] . 

Example: Here is an example of a discrete finite Gabor frame that is a nearly 
optimal Grassmannian frame in C m (see [2] for details about Gabor frames). 
Let m be a prime number > 5 and set g(n) = e 2mn l m for n = 0, . . . , m — 1. 
Then the frame {g^i}™^! where 

g k>l (n) = gin - k)e^ iln ' m , k, I = 0, . . . , m - 1, (25) 

satisfies \{gk,i, 9k',v)\ {0, l/^/m} for all gk,i ^ gk',i'i which follows from The- 
orem 2 in [I]. Hence M({g k ,i}Tr=o) = while (H2J yields 1/Vm + T as 
theoretically optimal value. Note that we can add the standard orthonormal 
basis to the frame {gkj}™^ without changing the maximal frame correlation 
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3 Infinite-dimensional Grassmannian frames 



In this section we extend the concept of Grassmannian frames to frames 
{fk}kLi m separable infinite-dimensional Hilbert spaces. As already pointed 
out in Section 121 the maximal correlation \(fk,fi) \ of the frame elements will 
depend crucially on the redundancy of the frame. While it is clear how to 
define redundancy for finite frames, it is less obvious for infinite dimensional 
frames. 

The following appealing definition is due to Radu Balan and Zeph Lan- 
dau [2J. 

Definition 3.1 Let {fk}kLi be a frame for 7i. The redundancy p of the 
frame {fk}kLi ^ s defined as 

•provided that the limit exists. 

Using the concept of ultrafilters Balan and Landau have derived a more 
general definition of redundancy of frames, which coincides of course with 
the definition above whenever the limit in (}2T)j) exists [2]. In this paper we 
will restrict ourselves to the definition of redundancy as stated in (J26|) since 
it is sufficiently general for our purposes. 

Remark: We briefly verify that the definition of frame redundancy by Balan 
and Landau coincides with our usual understanding of redundancy in some 
important special cases: 

(i) Let {fk}k=i be a finite frame for an m-dimensional Hilbert space 7i m . 
Let P : 7i n — ► 7i m denote the associated projection matrix with entries 
Pk,i = (fkiS-t-fi). We compute 

\JM ^— ^ / trace(.r) rank(_r) m 

which coincides with the usual definition of redundancy in finite dimensions. 

(ii) Let {g mt n}m,nez, where g mjn ( x ) — g (x — ma)e~ 2mnbx be a Gabor frame 
for X 2 (M) with time- and frequency-shift parameters a,b > 0. We have 
from [27] that 

(g m ,n, S~ l g rn:n ) = (g, S^g) = ab, for all m, n 6 Z, (28) 
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hence p = 1/ (ab) as expected. 

(iii) Assume {fk}kei is a uniform tight frame. Then (fk,S~ 1 fk) = l/A 
and therefore p = A, which agrees with the intuitive expectation that for uni- 
form tight frames the frame bound measures the redundancy of the frame ^U] ■ 

We need two more definitions before we can introduce the concept of 
Grassmannian frames in infinite dimensions. 

Definition 3.2 (|9|) A unitary system U is a countable set of unitary oper- 
ators containing the identity operator and acting on a separable Hilbert space 



Definition 3.3 LetlA be a unitary system and $ be a class of functions with 
1 1 V 1 1 2 = 1 f or V 9 e ^- We denote by J^iHyU, $) the family of frames {<fik}kei 
for TC of fixed redundancy p, such that 



We say that {fk}kex is a Grassmannian frame with respect to T{H,U, $) if 
it is the solution of 3 



for given p. 

In the definition above we have deliberately chosen $ such that it does 
not necessarily have to coincide with all functions in L 2 (J-C). The reason is 
that in many applications one is interested in designing frames using only a 
specific class of functions. 

In finite dimensions we derived conditions under which Grassmannian 
frames are also uniform tight frames. Such a nice and simple relationship 
does not exist in infinite dimensions. However in many cases it is possible to 
construct a uniform tight frame whose maximal frame correlation is close to 
that of a Grassmannian frame as we will see in the next theorem. 

The following definition is due to Gabardo and Han [TTj . 

3 For a frame {<fik}kei there always exists max{\ (tpk,<pi}\}, otherwise the upper frame 
bound could not be finite. 



n. 



Vk = U k (fo 



(p e $, U k eU, k el. 



(29) 




(30) 
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Definition 3.4 Let T denote the circle group. A unitary system U is called 
group-like if 

group(W) C TU := {tU : t E T, U E U}, (31) 

and if different U,V E U are always linearly independent, where group (U) 
denotes the group generated byU. 

Theorem 3.5 Let J-(Tl,L{,$>) be given, where U is a group-like unitary sys- 
tem. For given redundancy p assume that {fk}k&i is a Grassmannian frame 
for T(Tt,U,$) with frame bounds A,B. Then there exists a uniform tight 
frame {h k } k£ j with hk = Ukho, U k E U, such that 

jgKjtfc, ft,) | < JiS g (A , /,)| + 2max{|l - |1 - (32) 

Proof: Let S be the frame operator associated with the Grassmannian 
frame {f k }k£i- We define the tight frame {h k } k£ j via h k := ^fpS^^ f k . Since 
IA is a group-like unitary system it follows from (|29|) above and Theorem 1.2 
in ESI that 

(33) 
(34) 

(35) 



and 

frfc = VP S ~^ 2 fk = Vp~S~^U k f = y/pU k S~^f . 

Using Definition 13. II and (jUnj), we get (f k ,S~ 1 f k ) = - and therefore 



\hi 



p(S- l if k ,S- 1 if k )=p(f k ,S- 1 f k ) 



Vk El. 



Hence {h k } k <zj is a uniform tight frame. 
We compute 

KA,/i>|-IK/n>l| < K/*,/i>- I (36) 

< KfkJi-hN + Kh-h^h^l (37) 

< M + \\hi\\\\f k -h k \\, (38) 

where we have used the triangle inequality and the Cauchy-Schwarz inequal- 
ity. Note that 

\\fi-hi\\ = Wfi-VpS-^M (39) 

< W - VpS-i)\\\\M (40) 



< max < 1 1 



II 



(41) 
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Hence 



\{fhJi)\-\{h k ,hi)\ 



< 2 max < 1 1 — 




(42) 




and therefore 



max \(hk,hi)\< max \(f k , fi)\ + 2max<| |1 



Remark: (i) Although the canonical tight frame function h do not have to 
belong to $, it is "as close as possible" to the function / G <3>. Indeed, under 
the assumptions of Theorem 13.51 the (scaled) canonical tight frame {hk}k£i 
generated by h = ^S~2f minimizes ||/o — tpo\\ among all tight frames 
{tPk}kei hi U, L 2 (Ti.)) (in fact among all possible tight frames), cf. [23] 
and for the case of Gabor frames [2H]- However it is in general not true 
that {/ifcjfcgj also minimizes the maximal frame correlation max^j |(<^fc, ¥>z)| 
among all tight frames {(fik}kei- 

(ii) If the Grassmannian frame {fk}k^i is already tight, then the frame 
bounds satisfy A = B = p and the second term in the right-hand-side of (}3"2*j) 
vanished, as expected. 

(iii) Frames that satisfy the assumptions of Theorem l3.5l include shift-invariant 
frames, Gabor frames and so-called geometrically uniform frames (see 
for the latter). 



3.1 An example: Grassmannian Gabor frames 

In this section we derive an example for Grassmannian frames in X 2 (IR). We 
consider Gabor frames in L (R) generated by general lattices. 

Before we proceed we need some preparation. For x, y G M. we define the 
unitary operators of translation and modulation by T x f(t) = f(t — x), and 
M w f(t) = e 2mu,t f '(t) , respectively. Given a function / G L 2 (R) we denote the 
time-frequency shifted function f XtU by 

f« u {t) = e Mut f{t-x). (44) 

A lattice A of R 2 is a discrete subgroup with compact quotient. Any 
lattice is determined by its (non-unique) generator matrix L G GL(2,M) via 
A = LI? . The volume of the lattice A is vol (A) = det(L). 
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For a function (window) g E L 2 (R) and a lattice A in the time-frequency 
plane IR 2 we define the corresponding Gabor system Q(g, A) by 



g(g,A) = {M u) T x g, (x, u) E A} 



(45) 



Setting A = (x,u>) we denote g\ = MJT x g. If Q(g, A) is a frame for L 2 (IR) 
we call it a Gabor frame. As in remark (ii) below D efinit ion 13 . 1 1 we conclude 
that the redundancy of Q(g, A) is p = l/vol(A). A necessary but by no means 
sufficient condition for Q(g, A) to be a frame is vol(A) < 1, cf. [221- It is clear 
that max A ^v \ (g\,9x')\ will depend on the volume of the lattice, i.e., on the 
redundancy of the frame. The smaller vol(A) the larger max^y \ {g\,g\>)\- 

One of the main purposes of Gabor frames is to analyze the time-frequency 
behavior of functions ^3] . To that end one employs windows g that are well- 

1 2 

localized in time and frequency. The Gaussian f a {x) = (2/a)ie~ 7TCTX ,a > 0, 
is optimally localized in the sense that it minimizes the Heisenberg Uncer- 
tainty Principle. Therefore Gabor frames using Gaussian windows are of 
major importance in theory and applications. Our goal is to construct Grass- 
mannian Gabor frame generated by Gaussians. Recall that Q(f a , A) is a Ga- 
bor frame for X 2 (1R) whenever vol(A) < 1, see (221- Thus in the notation of 
Definition 1531 we consider H = L 2 (R), U = {T x M y ,x,y E A with vol(A) = 

12/ 

p}, and $ = {(f a | <Pa(x) = (2a)ie~ 7TX /<T , a > 0}. That means for fixed re- 
dundancy p we want to find A° among all lattices A with vol (A) = p and (p° 
among all Gaussians (p a such that 



is minimized. 

Since (p a = ipy a we can restrict our analysis to Gaussians with a = 1, as 
all other cases can be obtained by a proper dilation of the lattice. To simplify 
notation we write if := f\. 

Since T x and M u are unitary operators there holds | (g\, gy) \ = \ (g, gx-\) \ 
for any g E L 2 (K). Furthermore 1(^,^)1 is monotonically decreasing with 
increasing ||A|| (where ||A|| = a/|x| 2 + \oo\ 2 ) due to the unimodality, symme- 
try, and Fourier invariance of <p. These observations imply that our problem 
reduces to finding the lattice A° of redundancy p such that max | ((p, fx) I is 
minimized where 



max|((</? CT )A,(</v)A')l 



(46) 



A E {Le x , Le 2 , with e x = [1, 0] T , e 2 = [0, 1] T }. 



(47) 
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The ambiguity function of / G -L 2 (M) is defined as 



-t-oo 

Af(t, lj) = J f(x + t)f(x - \ )e~ 2 ^ dx. (48) 

— oo 

There holds \ (f,g)\ 2 = \ (Af,Ag) \ for f,g G L 2 (R). It follows from Proposi- 
tion 4.76 in that Aip is rotation- invariant. Furthermore Atp\ is rotation- 
invariant with respect to its "center" A = (x,u) which follows from (4.7) 
and (4.20) in j22j and the rotation-invariance of Aip. 

Next we need a result from sphere packing theory. Recall that in the 
classical sphere packing problem in M. d one tries to find the lattice A° among 
all lattices A in M. d that solves 

f Volume of a sphere 1 . , 

max { TTTT^ f ■ 49 

a [ vol(A) J v ' 

For a given lattice A the radius r of such a sphere is 

r = \{ min {||A-A'||}). (50) 
2 A, A e A 
A ^ A' 

Hence solving (j4^|l is equivalent to solving 

max{ min {||A — A'||}} subject to vol(A) = p (51) 
A^ A 

for some arbitrary, but fixed p > 0. Obviously the minimum has to be taken 
only over adjacent lattice points. 

Due to the rotation-invariance of Aip and A(p\ and since Aip(x, cu) is 
monotonically decreasing with increasing (x, uj) we see that the solution of 

min max (Kv^V^a)!} subject to vol(A) = p (52) 

A A as in 1471 

is identical to the solution of (|51|). It is well-known that the sphere packing 
problem (|5*T|) in R 2 is solved by the hexagonal lattice A hex , see jH]. Thus 
for given redundancy p > 1 the Gabor frame G(<p, Ah ex ) is a Grassmannian 
frame, where the generator matrix of Ah ex is given by 



J hex 



V2 1 

o 31 



(53) 
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Remark: (i) The result can be generalized to higher dimensions, since the 
ambiguity function of a <i-dimensional Gaussian is also rotation-invariant. 
Hence a Grassmannian Gabor frame with Gaussian window is always asso- 
ciated with the optimal lattice sphere packing in M? d . However in higher 
dimensions explicit solutions to the sphere packing problem are in general 
not known. jS]. 

(ii) If we define the Gaussian with complex exponent o = u + iv with u > 
(i.e., chirped Gaussians in engineering terminology) then it is not hard to 
show that a properly chirped Gaussian associated with a rectangular lattice 
yields also a Grassmannian Gabor frame. 

4 Applications of Grassmannian frames 

In this section we describe two applications of Grassmannian frames. The 
first one concerns wireless communication and involves infinite-dimensional 
Grassmannian frames, while the second one concerns coding theory and in- 
volves mainly finite Grassmannian frames. Another important application of 
finite Grassmannian frames is described in [21] where this concept is used to 
construct spreading sequences for CDMA. 

4.1 Grassmannian Gabor frames and OFDM 

We briefly describe an application of Grassmannian Gabor frames in wireless 
communication. Orthogonal frequency division multiplexing (OFDM) has 
emerged as attractive candidate for 4-th generation wireless communication 
systems. We refer to [15j for details about OFDM and to 39j for details 
about the connection to Gabor theory. The transmission functions of an 
OFDM system are 

(p kjL (t) = <p(t - kT)e 2nilF \ fcGZ,/ = 0,...,M-l, (54) 

where <p is a given pulse shape. To minimize distortions caused by additive 
white Gaussian noise the functions (fk,i are usually chosen to be mutually or- 
thogonal. OFDM transmission functions form a critically sampled or under- 
sampled Gabor system, since perfect reconstruction at the receiver requires 
TF > 1. The duality relations for Gabor systems provides the connection 
between Gabor frames and undersampled Gabor systems |14j . 
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The mobile wireless channel H is time-dispersive as well as frequency- 
dispersive which can lead to intersymbol interference (ISI) and interchannel 
interference (ICI). When the pulses (pk,i pass through the channel H, the 
mutual orthogonality is lost. Thus to mitigate ISI and ICI and to allow for 
a simple receiver structure it is important that the pulse shapes <fk,i are well 
concentrated in time and frequency. Since the Balian-Low theorem prohibits 
well- localized Gabor systems when TF = 1 one usually chooses TF > 1. In 
this situation Gaussians are a natural starting point for pulse shape design. 

If the pulses are well-localized in time and frequency we have 

\(Hip^i,Hip k±lt i ±1 )\ > \{H<Pk,i,tliPk±k',i±i')\, for k',l' > 1, (55) 

hence it is sufficient to consider the reduction of interference between pulses 
that are adjacent in the time-frequency domain. 

We can of course always increase the parameters T and F in order to re- 
duce interference, however this results in a lower data rate. Therefore in [U3J 
the OFDM scheme has been generalized by replacing the rectangular time- 
frequency lattice {(kT, lF)}k,iez by a general lattice, yielding a Gabor system 
of the form (jJ3J). As follows from Subsection 13. II when using Gaussian pulses 
in combination with a hexagonal-type lattice we can minimize the interfer- 
ence between pulses by maximizing the distance between adjacent lattice 
points without reducing the data rate. Gaussian pulses are not mutually or- 
thogonal, but we can easily orthogonalize this system via the "inverse square 
root method", see The "orthogonalized" pulses will have a somewhat 
larger interference as compared to Gaussian pulses, but this increase in in- 
terference can be estimated conveniently via Theorem 13.51 We refer to 40 
for more details. 

4.2 Erasures, coding, and Grassmannian frames 

Recently finite frames have been proposed for multiple description coding for 
erasure channels, see (2111201 El- We consider the following setup. Let {fk}k=i 
be a frame in E m . As in (jHJ) and (J7J) we denote the associated analysis and 
synthesis operator by T and T* respectively. Let / G E m represent the data 
to be transmitted. We compute y = Tf e E^ and send y over the erasure 
channel. We denote the index set that corresponds to the erased coefficients 
by 8 and the surviving coefficients are indexed by the set 1Z. Furthermore 



20 



we define the N x N erasure matrix Q via 

(0 ifk^l, 

Qki = I if k = I and k G £ , (56) 
[ 1 if = Z and k &1Z. 

Let £ represent additive white Gaussian noise (AWGN) with zero mean and 
power spectral density a 2 . The data vector arriving at the receiver can be 
written as y := Qy + e. 

The frame {fk\k=i is robust against e erasures, if {fk\k&n is still a frame 
for E m for any index set 1Z C {0, . . . , N — 1} with \1Z\ > N — e. In this case 
standard linear algebra implies that / can be exactly reconstructed from y 
in the absence of noise 4 . 

In general, when we employ a minimum mean squared (MMSE) receiver 
we compute the (soft) estimate 

f = (T^T n + a 2 I rn r l T^y, (57) 

where T-n is the analysis operator of the frame {fk}ken- This involves the 
inversion of a possibly large matrix (no matter if noise is present or not) that 
can differ from one transmission to the next one. The costs of an MMSE 
receiver may be prohibitive in time-critical applications. Therefore one often 
resorts to a matched filter receiver which computes the estimate 

f = cT*y, (58) 

where T is the analysis operator of the original frame {fk\k=i an d c > is 
a scaling constant. The advantage of an MMSE receiver is the better error 
performance while a matched filter receiver can be implemented at lower 
computational cost. 

Robustness against the maximal number of erasures is not the only per- 
formance criterion when designing frames for coding. Since any transmission 
channels is subject to AWGN, it is important that the noise does not get 
amplified during the transmission process. Yet another criterion is ease of 
implementation of the receiver. It is therefore natural to assume {fk\k=i 

4 For instance so-called harmonic frames are robust against up to N — m erasures |21| . 
which does not come as a surprise to those researchers who are familiar with Reed-Solomon 
codes or with the fundamental theorem of algebra. 
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to be a uniform tight frame, since in case of no erasures (i) the MMSE re- 
ceiver coincides with the matched filter receiver and (ii) AWGN does not get 
amplified during transmission. 

Our goal in this section is to design a uniform tight frame such that the 
performance of a matched filter receiver is maximized in presence of erasure 
channels. In other words the approximation error ||/ — /|| is minimized, where 
/ is computed via a matched filter receiver, i.e., / = jjT*y with y = Qy + e. 

We estimate the reconstruction error via 

~ Tfi 

\\f-f\\=\\f-^T*(QTf + e)\\ (59) 

777 777 

<-\\T*Tf - T*QTf\\ + -\\T*e\\ (60) 

777 

<^l|rwy/ll + *, (6i) 

where we used the notation P = I — Q. Since P is an orthogonal projection 
and therefore satisfies PP* = P there holds 

\\T*PT\\ 2 = \\(PT)*PT\\ 2 = \\PT(PT)*\\ 2 = \\PTT*P\\ 2 . (62) 

Hence we should design our tight uniform frame {fk}kex such that ||PTT*P|| 2 
is minimized, where the minimum is taken over all matrices P = I — Q, with 
Q as defined in 

Recall that TT* = {{f h f k )}k,iex, hence PTT*P = {{fi, fk)}k,iex r Fur- 
t her more 

\\PTT*P\\ 2 < > /||P7T*P|| 0O ||PIT*P|| 1 = maxy;K/ fc ,/ I )|. (63) 

lex. 

This suggests to look for frames for which max \(fk, fi)\ is minimized, in 

k,l,k^l 

other words to look for optimal Grassmannian frames. 



Remark: (i) In case of one erasure it has been shown in [21] (in case of un- 
known a) that uniform tight frames are optimal with respect to minimizing 
the influence of AWGN when using the MMSE receiver. In case of one era- 
sure uniform tight frames also minimize the reconstruction error when using 
a matched filter receiver. 

(ii) Holmes and Paulsen have shown that Grassmannian frames are optimal 
with respect to up to two erasures j2H]- This can be easily seen by minimizing 
the operator norm of the matrix PTT*P, which in this case reduces exactly 
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to the problem of minimizing max \ (fk, fi)\ for all k, I with k ^ I. 
(hi) There is strong numerical evidence that the optimal Grassmannian 
frames of part (b) in Corollary 12 . 61 are even optimal for three erasures (how- 
ever this is not the case for the frames constructed in part (a)), 
(iv) Grassmannian frames are in general not robust against N — m erasures 
if iV > m + 2. 

Example: We elaborate further an example given in [2T], where the authors 
consider the design of multiple description coding frames {fk}k=i m ^ m with 
m = 3 and N = 7. As in Examples 4.2 and 4.3 in [21] we consider an 
erasure channel with AWGN, but unknown noise level. Without knowledge 
of a the the reconstruction formula of the MMSE receiver simplifies to / = 
(T^T^) -1 !^?/. Standard numerical analysis tells us that the smaller the 
condition number of T^T-ji the smaller the amplification of the noise in the 
reconstruction [T^]. We therefore compare the condition number of different 
uniform tight frames for m = 3, N = 7 after up to four frame elements have 
been randomly removed. 

We consider three types of uniform tight frames. The first frame is an 
optimal complex-valued Grassmannian frame, the second one is constructed 
by taking the first three rows of a 7 x 7 DFT matrix and using the columns of 
the resulting (normalized) 3x7 matrix as frame elements (this is also called 
a harmonic frame in [21])- The last frame is a randomly generated uniform 
tight frame. Since all three frames are uniform tight, they show identical 
performance for one random erasure, the condition number of the frame 
operator in this case is constant 1.322. Since the frames are of small size, 
we can easily compute the condition number for all possible combinations 
of two, three, and four erasures. We then calculate the maximal and mean 
average condition number for each frame. As can be seen from the results 
in Table 14.21 the optimal Grassmannian frame outperforms the other two 
frames in all cases, except for the average condition number for two erasures, 
where its condition number is slightly larger. This example demonstrates the 
potential of Grassmannian frames for multiple description coding. 
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